clc;
clear;
f=@(x)cos(15*x);
tic,
I1=integral(f,0,1000,'RelTol',1e-20),
toc
disp('-----------------------------')
syms x;
tic,
I=int(cos(15*x),x,0,1000),
vpa(I)